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Abstract 

The CELESTE atmospheric Cherenkov detector ran until June 2004. It has ob- 
served the blazars Mrk 421, lES 1426+428 and Mrk 501. We significantly improved 
our understanding of the atmosphere using a LIDAR, and of the optical through- 
put of the detector using stellar photometry. The new data analysis provides better 
background rejection. We present our light curve for Mrk 421 for the 2002-2004 sea- 
son and a comparison with X-ray data and the 2004 observation of lES 1426+428. 
The new analysis will allow a more sensitive search for a signal from Mrk 501. 
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1 Introduction 



CELESTE was a Cherenkov experiment using 53 heliostats of the former 
Electricite de France solar plant in the French Pyrenees at the Themis site. 
It detected Cherenkov light from electromagnetic showers produced in the at- 
mosphere by the 7-rays coming from high energy astrophysical sources. The 
light is reflected to secondary optics and photomultipliers installed at the top 
of the tower. Finally it is sampled to be analysed (Pare 2002). 

To constrain the energy scale of the experiment we improved the optics simu- 
lation, now in good agreement with the data. The data analysis has also been 
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improved so that we have better background rejection. We present the hght 
curve for Mrk 421. 



2 Constraining the energy scale 

The simulation has been reexamined to reduce the uncertainties on the energy 
scale of the experiment (Brion 2003). The LIDAR operating on the site for 
atmospheric monitoring provided a better determination of the atmospheric 
extinction (Bussons Gordo 2004). A stellar photometry study, focussing on 
the comparison between simulations and data on bright stars' currents, has 
been done on the star 51 UMa (Mb = 6.16) which is in the field of view 
(FOV) of Mrk 421. This showed that the old simulation was too optimistic. 
All mirror reflectivities were decreased subsequent to new measurements. The 
nominal focussing of the heliostats was degraded after a study of star image 
sizes. We also verified the photomultiplier gains. The results for the old data 
set (40 heliostats, new is 53 helisotats, see § 3) are presented in figure 1: the 
effect of these changes is smaller for 7-ray showers (an extended light source) 
than for stars (point sources). 
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Fig. 1. ON-OFF illumination from the star 51 UMa (Mb = 6.16) in the FOV of 
Mrk 421 as a function of hour angle for pointing at 11 km: the new simulation 
with our corrections (red stars) fits the data well (black squares) whereas the old 
simulation (blue circles) was 50 % too high. 

This study helped us to perform our selection criteria for the data: all data 
that are too low in currents on the star have also low trigger rates. In order 
to trigger CELESTE, the heliostats, except the veto (see § 3.1), are split into 
6 groups. For each of them, the analog sum of signals gives the first level 
trigger. Then, a logical pattern is defined on the majority of the triggering 
groups. Thus, in the case of the source Mrk 421, we reject all data with trigger 
rates under 20 Hz for the old data set and under 16 Hz for the new data set. 

We also looked at the proton rates as a standard candle for the detector 
which should be stable for good quality nights. These rates are determined 
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with high offline threshold cuts to avoid trigger bias, and are therefore low 
(typical trigger rates ~ 22 Hz). We've shown that they are correlated with 
the currents on star 51 UMa for selected data (new data set, figure 2 (a)). 
The data with low rates are now still rejected but perhaps some of them could 
be corrected. Indeed, three doubtful zones. A, B and C on the figure, can be 
distinguished. The zone A can be interpreted as bad nights with thick cloud 
cover (weak star, little Cherenkov light), the zone B as nights with aerosols 
and above average extinction of Cherenkov and starlight, and the zone C as 
nights with high clouds that stop starlight but not Cherenkov light. The data 
in this last zone may therefore be used. Figure 2 (b) shows the correlation 
between the proton rate and the trigger rate for the same data set. Defining 
a selection criteria for each type of acquisition, and based only on these rates, 
would be very interesting for sources that don't have any star nearby for a 
photometry study. 
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Fig. 2. (a) Proton rate as a function of ON— OFF current for heliostat EOS that 
sees 51 UMa when pointing Mrk 421. For a typical photomultiplier tube gain of 
5.6 X 10'*, 2.5 |xA corresponds to 0.28 p.e./ns, as seen in figure 1. (b) Proton rate as 
a function of trigger rate for data on Mrk 421. 



3 Analysis improvement 



Since the 2000 status of the experiment (de Naurois 2002) three main changes 
were made for the analysis. First, the selection criteria of the data are stricter 
for the current and trigger rate stability, and as we've shown for the trigger 
and proton rate value (work in progress). The experiment has been upgraded 
from 40 to 53 heliostats. We use part of them to broaden our narrow FOV. 
Finally, we have found a new method to exploit the FADC information to 
reject the hadronic background (Manseri 2004 a, b). The second and third 
point are developed hereafter. 



3 



3.1 The veto configuration 



Hadronic showers have a more chaotic and extended development than electro- 
magnetic showers. To measure the extent of the shower, we artificially broaden 
the FOV: as before, all heliostats aim at 11 km above the ground in the direc- 
tion of the source, where the maximum of the shower is supposed to occur in 
our energy range. But 12 heliostats, distributed around the edge of the field, 
sample a ring of 150 m around that point (figure 3). 

Because of the compactness of electromagnetic showers, the light does not 
illuminate these 12 heliostats, named veto, contrary to hadronic showers (fi- 
gure 4). So we require that no veto be illuminated. 




Fig. 3. Scheme of the veto configura- Fig. 4. Distribution of the number 

tion. of illuminated veto heliostats (simu- 

lations). 



3.2 FADC information 



To detect low energy 7-rays, we use the sum of the individual digitized signals 
to increase the signal-to-background ratio. The summation includes a correc- 
tion for the sphericity of the wavefront, assumed to be centered in the 11 km 
plane. Assuming a wrong position for this center (impact parameter) broadens 
the sum: the height-over- width ratio, (H/W), decreases. We compute (H/W) 
for different assumed positions. The impact parameter is the position for which 
the (H/W) ratio is maximum, denoted by {H/W)max- This is valid for 7-rays 
(figure 5 (a)) but not for protons for which the wavefront is not spherical (fig- 
ure 5 (b)). A measurement of the flatness of these 2D-distributions is given 
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by the following estimator: 



^^^^^^^ [jmu-., 



over 24 positions 



where {H /W)2QQm is the average of {H/W) over 24 positions along a ring 
200 m from the maximum position. For 7-rays there is a clear maximum and 
this estimator takes low values. For hadrons, it is usually larger (figure 5 (c)). 
We have applied ^ < 0.35. 
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Fig. 5. Left: the shape of the FADC sum depends on the position of the centre of the 
wavefront we assume. Right: views of the height-over- width ratio (H/W) computed 
for discrete positions in the 11 km plane (a) for a 100 GeV 7 ray and (b) for a 500 
CeV proton, (c) Distribution of the ^ parameter for 7-rays and OFF data. 

A cut on this estimator is the biggest contribution to an improvement in Crab 
sensitivity from 2.2 cr/-\/h to 5.8 a/\/h.. The Crab detection is also stable 
(figure 6). 



4 Blazar observations 



This stable analysis and good sensitivity provide a good detection of the 
Mrk 421 fiares. A 19 a detection during 10 h since October 2002 gives a 
mean of 5.6 7/min. The light curve is presented in figure 7. With the new 
acceptances we will deduce a spectral measurement for Mrk 421 with smaller 
uncertainties. 

The source lES 1426+428 is still not detected with our new analysis during 
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Fig. 6. Stability of the Crab detection after hour angle efficiency corrections. Line 
is 6.5 =b 0.5 7/min. 
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Fig. 7. Light curve of Mrk 421. Left: since 2002 seen by CELESTE (red points in 
7/min corrected for hour angle efficiency, 1 point is 20 min data, Crab transit rate 
shown for reference) and RXTE/ASM (blue squares in mCrab, 6 day bins). Right: 
light curve of Mrk 421 zoomed on the 17*^ of March 2004 flare seen by CELESTE. 

the March 2004 observation for 4.4 h data. Since the data have been taken in 
the same conditions as for the Crab, we can use the sensitivity of 5.8 a/Vh to 
remark that for a 3 a observation for these 4.4 h data, the flux /iesi426 would 
have to be /iesi426 = 5 8^^4.4 f'^^^^ ~ Ic^_ \ better upper hmit determination 
is in preparation. 

Finally, the blazar Mrk 501 has been observed during 14.5 h with the pre- 
vious experiment at 40 heliostats. The old analysis gave a 2.5 a signiflcance 
(with 850 photons which would give 1 7/min). The new analysis will allow 
an improved investigation of Mrk 501 's behaviour in this unexplored energy 
range. 
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